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This  paper  presents  a  proof  of  an  interpolation  result  related  to  the  approx¬ 
imation  theory  for  higher  order  finite  element  or  spectral  methods  when  C*  (or 
higher)  regularity  is  convenient  for  the  finite  dimensional  subspaces.  This  can 
be  a  natural  choice  for  example  for  Stokes  problem,  the  biharmonic  problem  or 
higher  order  plate-  and  shell  models.  We  show  that  one  gets  the  same  interme¬ 
diate  spaces  whether  one  1)  interpolates  between  two  Sobolev  spaces  defined  on 
a  domain  with  nonsmooth  boundary  first  and  then  enforces  the  homogeneous 
boundary  conditions  afterwards  or  2)  interpolates  between  two  Sobolev  spaces 
where  the  homogeneous  boundary  conditions  are  enforced  throughout  the  inter¬ 
polation  proces. 


Keyc_words.  Interpolation,-  Peetre."  boundary  conditions,"  nonsmooth  domains," 


small  angle  elliptic  regularity. 
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1.  Introduction 

The  aim  of  this  note  is  to  prove  an  interpolation  result  for  domains  in  R2  with  finitely 
many  corners  and  otherwise  smooth  boundary.  We  consider  a  bounded  open  set  0  of 
R2,  whose  boundary  is  a  curvilinear  polygon  of  class  C°°  (see  [6]).  We  denote  each 
of  the  C°°  curves  which  constitute  the  boundary  by  T,-  for  some  j  ranging  from  1  to 
N.  The  curve  TJ+1  follows  T,-  according  to  the  positive  orientation,  on  each  connected 
component  of  T.  We  denote  by  Cj  the  vertex  which  is  the  end  point  of  Tj  and  by 
0£j  the  measure  of  the  angle  at  Cj  (toward  the  interior  of  0).  By  a  corner  we  mean 
a  vertex  Cj  with  an  angle  oij  not  in  the  set  {0,7r,27r}.  The  result  is  an  extension  to 
R£  of  the  one  in  [3],  [1]  for  Hi  which  would  be  useful  in  approximation  theory  for 
Sobolev  spaces,  see  [7],  Remark  2.2.9,  [14],  Remark  4.2,  [2],  and  [11],  the  line  following 
(III. 26)  in  the  proof  of  Thm.  III. 2.  For  example,  consider  solving  Stokes  problem  via 
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the  p  version  of  the  finite  element  method  or  a  polynomial  spectral  method.  Then  the 
discrete  velocity  Up  is  an  elliptic  projection  onto  a  finite  dimensional  subspace  Zp  of 
Z  =  [ //q ( Q ) ] 2  D  Ker(div)  centering  interest  on  the  approximation  problem.  Introducing 
stream  functions  ( U  =  rot <j>,Up  =  rot^p)  will  translate  this  approximation  problem  to 
Hq(Q,).  Now  one  gets  for  free  an  energy  estimate  -  \\f>  -  <^p||2  bounded  when  <f>  £  Hi 
only  -  and  there  exists  constructive  approximation  estimates  -  \\<f>  —  (f>p || 2  <  Cp'-'Mt 
when  <j>  €  fl  Hl(£l)  for  t  >  7/2,  see  [14].  Now,  one  wishes  to  interpolate  between 

these  spaces  and  hopes  to  get  spaces  that  coincide  in  some  sense  with  the  ones  predicted 
by  regularity  theory,  but  the  trace  constraints  on  a  nonsmooth  boundary  makes  this 
identification  nontrivial.  In  general  such  an  identification  is  useful  for  higher  order 
finite  element  or  spectral  methods  when  C 1  (or  higher)  regularity  is  convenient  for  the 
finite  dimensional  subspaces.  This  can  be  a  natural  choice  for  example  also  for  the 
biharmonic  problem  or  higher  order  plate-  and  shell  models. 

Let  H‘(Q)  be  the  standard  Sobolev  space  of  order  s  based  on  L2  with  corresponding 
norm  ||  •  ||,.  is  the  set  of  functions  in  i/m(fl)  for  which  the  traces  of  the  function 

and  its  normal  derivatives  up  to  order  m  —  1  vanish  on  dfl. 

We  shall  use  the  interpolation  spaces  of  Peetre  (see  e.g.  [4])  in  the  cases  1  <  q  <  00 
where  we  define  [#‘(0)  fl  explicitely:  For  u  6  Hl  fl  H™,  we  set 

K(u,s)  =  inf  (IML  +  s|M|t) 

u  —  v  4-  to 

v  €  H™,  w  €  Hl  D  Hf 


and  we  define  the  norm 


u 


\9,q 


Then 


[ff*(fl)  n  H?(n)]e,q  =  {ue  isrnn)  =  IMIlv]#>,  <  00} 

[Hl(U),  Hm( fl)]®,,  is  defined  similarly.  Note  that  this  space  will  be  a  Sobolev  space  if 
we  choose  q  =  2  and  in  general  a  Besov  space. 

2.  The  interpolation  result 

We  state  and  prove: 


ed  □ 

Proposition  1  Let  fl  C  R2  be  piecewise  C°°  with  finitely  many  comers  of  angles  in  tion - 


(0,  27t)\  {7r}.  Then  the  following  identity  holds  for  all  6  6  (0, 1),  1  <  q  <  oo  and  t  >  m, 
t  i  m  +  {^,...,m  -  i},  m  €  Z+: 

[//‘(fl)  n  f/0m(fi),  ff0m(n)]*,,  =  [jy*(n),  Hm(Q)}e,q  n  h?(H) 
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Proof:  We  follow  the  main  ideas  of  [31  but  have  weights  be  unity  for  simplicity,  see 
also  [15]  and  [1], 

The  inclusion  from  left  to  right  follows  directly  from  the  definition. 

The  proof  of  the  reverse  inclusion  can  through  a  partition  of  unity  be  reduced  to 
considering  a  domain  12  with  one  corner  of  angle  a  £  (0, 2ir)  \  {7r} .  We  shall  distinguish 
between  two  cases:  whether  a  £  (0,7r)  or  (7r,27r). 

Case  1:  a  £  (0,7r).  Then  there  exists  a  linear  transformation  from  12  to  12  with  a 
corner  of  angle  a  <  min{o>0,7r}  where  u0  will  be  introduced  in  the  next  section  as  a 
sufficiently  small  angle  that  a  certain  shift  theorem  will  hold.  Let  L  be  the  associated 
map  of  functions  defined  on  12  to  functions  defined  on  12.  If  w  £  #‘(12),  then  we  let 
v  =  PfiUi  denote  the  solution  to 

+  v  =  (  — A)mu>  -f  u;  in  12  .  . 

v  €  #£*(  12)  1 ’ 

Thus  Pp  is  a  projection  from  #m(12)  to  #£*(12).  As  proven  in  the  next  section  on 
regularity,  there  exists  u>0,  dependent  on  m  and  t,  such  that  the  following  shift  theorem 
holds  provided  a  <  lu0: 

p||jf‘  <  C||(-A)mu;  + 

In  particular,  Pn  =  L  o  P~  o  L~l  €  B(Hm(Q),  #£*(12))  and  Pn  6  £(#‘(12),  P‘(12)  D 
#£*(  12))-  Thus,  by  interpolation, 

Pn  G  £([#'(12),  #TO(  12)]s,„  [#‘(12)  n  #£*(12),  #£*(12)], „) 

Since  Pn|j/™(n)  =  I  (the  identity), 

[#'(«),  #m(fi)]«,,  n  #£*(12)  C  [#‘(12)  n  #£*(12),  #£*(12)],, ,  (2.2) 

Case  2:  a  E  (7t,27t).  Let  B  be  a  ball  centered  at  the  corner  and  containing  12. 
By  Pn,  we  denote  the  Stein  extension  ([13]  Chapter  VI,  Sec.  3)  of  functions  on  12  to 
functions  on  B  vanishing  at  dB.  Let  12e  =  B  \  12  and  let  Pp<  be  the  Stein  extension  of 
functions  on  12c  to  all  of  B.  Theorem  5  in  [13]  states  that  Pp  £  £(#fc(  12),  Hk(B))  and 
Pn<  <E  £(#*(  12c),  #*(£)),  V*  £  N.  Now  define 

Pn  =  #ne  °  Poc  0  #n  +  {I  —  #ne  °  £n) 

with  Pn«  being  the  the  same  operator  as  in  Case  1.  Then  Pn|jr0m(n)  =  I  and 

Pn  €  £((#‘(12),  #m(n)],,„  [#‘(12)  n  #£*(12),  #£*(12)],, ,) 

which  ends  the  proof  of  the  proposition. 

//// 

Remark  1:  We  have  explicitely  excluded  vertices  of  angles  0,7r,  or  2w.  In  these 
cases  it  is  not  possible  to  map  linearly  onto  a  domain  of  sufficiently  small  angle.  In 
case  12  were  a  polygon  the  exclusion  only  amounts  to  2ir. 

Remark  2:  The  theorem  and  proof  hold  for  conical  points  in  R3  almost  verbatim. 
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3.  Regularity  for  small  angles 

In  [lOj  it  is  stated  that,  given  k  £  N,  if  the  domain  contains  only  corners  of  sufficiently 
small  angles  and  /  £  H k,  then  the  solution  (u)  of  a  Dirichlet  problem  with  zero 
boundary  conditions  and  a  2 m  order,  elliptic  operator  ( Lu  =  /,  u  6  H ™)  belongs  to 
Hk+2m .  We  present  a  proof  here  following  and  expanding  upon  the  ideas  in  [10]  pp 
292-294  and  [5]  extending  to  the  case  where  /  £  Hk,  k  >  —  m.  We  use  the  notation  of 

[10]' 

Let  L  be  elliptic  of  order  2m  and  u  £  H™  be  the  solution  of 

Lu  =  f 

In  [10]  a  technique  is  used  that  involves  a  combination  of  1)  looking  at  L0( 0,  the 
principal  part  of  the  operator  L(x,  £x)  with  coefficients  fixed  at  the  origin,  2)  changing 
to  polar  variables  (r,  u/),  so  that  LqU  =  /  takes  the  form: 

4-'  7>*i 

«l.«3=0 

3)  making  the  change  of  the  radial  variable  (g  =  In  ~)  so  that  L0u  =  /  now  takes  the 
form 

2m 

£  = /•=-■”'  =  f 

1  de^dw* 

and  4)  taking  the  Fourier  transform  with  respect  to  the  “radial”  variable  (p).  The 
domain  then  consists  of  angles  u  £  D  -  an  interval  (for  R2  -  in  higher  dimensions  a 
cylinder  with  smooth  boundary).  The  final  form  of  L0u  =  f  is 

L0{u;,iX,--^-)u  =  F 

The  boundary  conditions  undergo  similar  transformations.  Let  /2(A)  be  the  resolvent 
operator  -  a  meromorphic  function  of  A  -  associated  with  the  resulting  boundary  value 
problem.  In  [10]  and  [9]  it  is  shown  that  if  /  £  //*,  then  u  £  Hm  has  the  expansion 


u  =  £  £a/'r  ’A,'log'r  £  P,-„(rlog9r)  (3.1) 

ho<ImAi<h,=0  «=° 

*0  Halloo 

+  £  clprllogpr  j>ip(u>)  +  w  (3.2) 

1= 0  p=0 

where  h0  —  —  1  +  m,  h  =  -1  +  A:  +  2m,  A j  are  the  poles  of  /2(A)  of  multiplicity  nj,  Pj,q 
are  polynomials  of  degree  [h  —  ImAj]  whose  coefficients  are  C°°  functions  of  uj  as  are 
iptp,  and  w  £  Hk+2m.  From  this  expansion  we  see  that  the  smoothness  of  u  depends  on 
the  poles  Aj  of  the  function  /2(A)  which  lie  above  the  straight  line  ImA  =  —  1  +  m.  We 
will  show  that, 
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Lemma  1  Given  any  positive  h,  there  exists  w0  such  that,  if  the  angle  of  the  comer  is 
smaller  than  u>0,  then  the  strip  —  1  +  m  <  ImX  <  h  contains  no  poles  of  R(X). 

Proof:  Let  h  be  given  and  A0  be  a  pole  of  R( A)  lying  in  the  strip  —  1  +  m  <  ImA  <  h. 
When  A  =  A0,  there  exists  a  nonzero  solution  u0(u>)  to 

Q 

Li(A,u;,  —  )u°  =  0  in  D 


u0  = 


du0 


& 


m- 1 


Uq 


=  0  on  r 


dtu  duim~x 

Now,  L\u0  =  Lqu0  +  XLiUo,  where  the  operator  L\  contains  derivatives  of  order  less 
than  2m.  Since  this  system  is  elliptic  for  all  real  A,  it  is  elliptic  for  A  =  0.  So  L0  is 
elliptic.  Let 

I(u )  =  J_(Liu)uduj 

which  at  u0  is  zero:  integrate  by  parts 

A(wo)  +  ^2(^0)  +  Isiy-  0)  = 

(  rr  t  sdmu0dmui 

Js{KJu)9^^ 


,  v~'  \*+i  t  \®m  *uo Ju0l 

+  E  A  ) 

0<i+j<2m  uw  UW 

+A2ma0o(w)u0:Uo}du; 

=  0 


By  ellipticity  of  L0, 

|Re/1(u0)|  >  a0  J-  \~~\2<L>  -  C0  J-  \u0\7du/ 

where  a0  and  Co  do  not  depend  on  u0  or  -  the  diameter  of  D. 

|/J(«o)l  <  e /-  |^l’^  + 

as  sketched  in  [10]  and 

|/3(«o)|  <  C|A|Jm  f  luofdw 
Jd 

where  C  does  not  depend  on  <r .  Also 
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Upon  contracting  like  terms,  substituting  this  last  inequality,  and  cancelling  \u0\2dw, 
we  get  from  fRe/^  <  |/2|  +  I/3I 

al(T-2m  <  1  4-  jAI2*"1"1)  +  |A|2m 


for  some  Oj  >  0  independent  of  a.  Since  ImA  £  (  —  1  +  m,h),  if  A  =  petB,  0  £  [— f ,  y), 
given  any  e  >  0,  there  exists  <r  sufficiently  small  so  that 


and  thus  2 md  £  (  — e,  e)  U  (2rmr  -  e,  2m7r  +  e)  =  (-e,  e)  such  that  Re(A2m)  >  0.  Having 

|ReJ*|  >  l\X\2m  L  \u0\2du  ~  <7(m)|IrnA|2m  f  \u0\2du> 

2  JD  JD 

If  we  again  contract,  substitute,  and  cancel  as  before,  but  now  in  iRe/x  +  Re/3|  <  |/2| 
(using  a0o  and  amm  have  same  signs  and  Re(A2m)  >  0),  we  get 

<r~2m  +  CJAI2"*  <  CalAI2^'1)  +  C3 


admitting  no  solutions  A  for  sufficiently  small  <r. 

nn 

Remark  3:  Another  way  of  proving  this  lemma  for  the  biharmonic  operator  is  by 
checking  that  one  can  choose  the  angle  a  such  that  the  roots  of  the  equation 

sinh2(ru;)  =  r2  sin2  u> 

see  [6]  (7, 2, 2,1),  except  for  -*  and  0  all  have  sufficiently  small  (negative,  large  absolute 
value)  imaginary  value.  Note  1  +  ir  =  — iA. 

Remark  4:  For  the  H2  interpolation  to  hold,  it  now  suffices  to  quote  [6],  Thm. 
7. 2. 2. 3.  The  interpolation  result  is  thus  essentially  -  along  with  the  reasoning  of 
the  proof  of  the  Prop,  said  a  localization  of  the  poles  of  R{ A)  for  the  biharmonic  - 
a  consequence  of  the  analysis  in  Grisvard’s  monograph  [6].  Such  an  analysis  was  first 
done  in  [8]  and  [12]  for  Stokes  problem  (which  via  the  stream  function  connects  to  the 
biharmonic  problem). 

For  m  >  2  we  shall  employ  a  recent  regularity  theorem  in  [5]. 

Lemma  2  Assume  that  k  >  —m,  k  $  —m  +  {5,  •  •  •  ,m  —  |}  and  let  f  £  Hh(U)  (ft 

is  defined  above).  Then,  for  sufficiently  small  a,  the  solution  u  £  to 

(—A)mu  +  u  =  f  in  ft 

belongs  to  Hh+2m(ft)  and  ||u||jf*+am  <  C\\f\\gk. 


(3.3) 
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Proof:  The  statement  is  really  a  corollary  of  Lemma  1  and  a  recent  shift  theorem  in 
[5],  Corollary  (5.2)  proven  in  [5]  §10.  It  is  stated  for  a  cone  in  Thm.  (1.11).  In  order 
to  apply  this  result  we  select  the  angle  a  sufficiently  small  that  Lemma  1  ensures  that 
iZ(A)  has  no  poles  in  the  strip  ImA  £  [m  —  l,t  —  2].  This  in  turn  implies  Dauge’s 
condition  (C2*)  (~  (R2))  as  follows:  If  — iA  £  N,  then  [5]  Cor.  (4.6’)  yields  (C2")  and 
if  —  i\  £  N,  then  Cor.  (4.9)  along  with  Cor.  (4.15)  and  the  fact  that  a  —  2  (see  [5]  p. 
39)  concludes  the  proof. 

till 

Remark  5:  For  m  >  2  it  is  possible  to  prove  Lemma  2  directly  from  [10]  when 
k  >  -1,  k  €  Z.  [10]  has  the  result  for  /  £  Hk  and  k  £  N.  Then  it  is  possible  to 
prove  for  /  £  Hk  when  k  >  —1  by  generalizing  the  trace  theorem  7. 2. 2. 3  in  [6]  for  a 
domain  Q  with  one  corner  C  of  sufficiently  small  angle  a  between  the  two  linear  pieces 
Tj,  j  =  1,2.  First,  by  this,  one  finds  v  £  Hk+2m(D)  0  H™(Q)  such  that 

(-A )mv  +  V  -  f  e  Hk{ti)  (3.4) 

Then  one  applies  to  w  =  u  —  v,  Kondrat’evs  result  (when  /  £  Hk)  with  Lemma  1 
choosing  a  sufficiently  small  that  w  £  Hk+2m  where  v  £  Hk+2m  is  the  solution  to  (3.4). 
In  [6]  the  generalization  of  Kondrat’ev’s  weighted  spaces  is  given  for  k  =  —  1.  It  seems 
difficult  however  to  go  to  all  the  remaining  negative  integers. 
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